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Abstract. Recently, Biswas and Milovic [Appl. Math. Comput. 208 (2009) 209-302] 
have found optical one-soliton solutions of a fourth order dispersive cubic-quintic non- 
linear Schrddinger equation. In this comment, we first show there are mistakes in 



GO 

<**] ■ the paper and demonstrate that the obtained solutions do not satisfy the considered 

equation. And then we reconstruct a series of analytical exact solutions by means of a 
direct ansatz method and F-expansion method. These solutions include solitary wave 

(N 

solutions of the bell shape, solitary wave solutions of the kink shape, and periodic wave 

m ' 

i solutions of Jacobian elliptic function. 

q 

^ ■ 1. Analysis of the solutions given in Ref. [1 
O 

As is well known, the investigation for solition solutions of nonlinear Schrodinger equation 
is always an important and attractive topic. Very recently, Biswas and MilovicpQ considered 
the higher order dispersive cubic-quintic nonlinear Schrodinger equation, 

iq t + aq xx - bq xxxx + c(|<?| 2 + d\q\ 4 )q = 0, (1) 

and obtained the optical soliton solution of Eq.(fT]). However, we find there are mistakes in 
the paper pQ and the obtained solution does not satisfy Eq.(TjQ). 

Biswas et al.pQ first introduced the transformation, 

q(x,t) = P{x,t)e i( - KX+ujt+e \ (2) 

where P(x,t) is a real function to be determined later, and k,uj are real constants. By using 
the transformation ([2]), Eq.([T]) is converted into a complex differential equation of P(x,t), in 
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which the real and imaginary parts read, 

dP dP d 3 P 

__ 2K((1 + 2 ^ ) _ +4tK - ?= 0, (3) 

(w + ok 2 + bK A )P - cP 3 - cdP 5 - (a + 6bK 2 ) + = 0. (4) 

ox ox 

Then the solution of Eqs.(J3])-(|l]) was supposed as 

p = 7\ r^J> r = B(x-vt). (5) 

(A + coshr) p ' y ' w 

Substituting Eq.© into Eqs.([3|)-(jl|), the authors obtained two expressions with respect to 
cosIit and sinhr. We have noticed that there are many mistakes about the expressions 
(12)-(13) given in Ref.p. Equating the coefficients of 1/(A + coshr) p+ -?(j = 0, • • • ,4) of the 
obtained expressions, the values of A, B,uj, A and v were determined. 

At last, the authors obtained the optical soliton solution of Eq.© as follows, 

ffOM) = 1 Tt4t ii-KM+9) (6) 

Hy ' ' A + cosh( J B(x -vt)) y 1 

where A, B,u), A and v were given by Eqs.(16)-(21) of Ref.[T]. 

However, the "solution" @ does not satisfy Eq.(TlJ. We can note this fact without substi- 
tuting © into Eq.(TjQ). The solitons are the results of a delicate balance between dispersion 
and nonlinearity, thus it is impossible that the linear partial differential equation (|3|) admits 
the bell type solitary wave ([5]). 

To be on the save side we have substituted Eq. © with p = 1 into Eq. ([3]) and have obtained 
the following expression, 



v + 2qk + 4bK 3 - 4bnB 2 24bnB 2 cosh r 246k5 2 sinh 2 r) 



+ 



A B sinh r. 



(A + cosh t) 2 (A + cosh r) 3 (A + cosh r) 4 

We can see that this expression is equal to zero only in two cases. One is A = or B = 0, 
and the other is k = v = 0. This means that the "solution" obtained by Biswas et al. in 
PQ is not correct. 

2. New optical solitary wave solution of Eq.fll]) 

In the following, we adopt the ansatz solution of Li et al.[2] in the form 

q(x,t) = E{x,t)e i{ - kx - wt+e \ (7) 



2 



where E(x,t) is the complex envelope function, and k, uj are real constants. Substituting 
Eq.® into Eq.([T|) and removing the exponential term, we can rewrite Eq.([T|) as 
iE t + 2ik{a + 2bk 2 )E x + (a + 6bk 2 ) E xx - AibkE xxx - bE xxxx 

(8) 

+ (uj - ak 2 -bk 4 )E +c\E\ 2 E + cd\E\ 4 E = 0. 

We now take the complex envelope ansatz function E(x,t) as 

E(x,t) = if3 + A tanh(f), £ = px + st, (9) 

where P, X,p, s are real constants. Substituting Eq.© into Eq.© and setting the coefficients 
of tanh(£) J (j = 0, • • • ,5) to zero, one obtains the following algebraic equations: 

A (dcX 4 - 24 p 4 b) = 0, A(A 3 /3 cd + 24 bp 3 k) = 0, 

A(s - 2 A p 3 cd - A c + 2 kpa + 32 &p 3 /c + 4 6p/c 3 ) = 0, 

A (A 2 c + 2 A 2 /3 2 c<i + 40p 4 b + 2p 2 a + 12 /c 2 p 2 &) = 0, 

A (2p 2 a + k 2 a - p 4 cd - f3 2 c + k 4 b + Wp 4 b - u + 12 A; 2 p 2 6) = 0, 

p bk 4 - u - P 3 c - 2 a P X k - p 5 cd - sA - 4 bpX k 3 - 8 bp 3 kX + /3 ak 2 = 0. 
Solving it we obtain one set of nontrivial solution, 

s = 8bpk(k 2 +p 2 ), = -M, w = 2p 2 a + 3/c 2 a + 37£; 4 & + 52/cV& + 16p 4 6, 

2p 2 (30 6fc 2 + 20 6p 2 + a) 12frjp 2 ^ 

C " A 2 ' A 2 (30 6/fc 2 + 20 6p 2 +a)' 

From © , © and (fTUj) , we obtain the optical solitary wave of Eq. ([1]) , 
ft ( Xjt ) = (_i|A + Atanh(px + 86pA;(fc 2 +p 2 )t)) e i(^-(2 P 2 a+3fc 2 a + 37^ +52fc Vfe+16 P ^) t+ ,) 



p 

where p, are determined by the last two identities of Eq. (|10p . From ([9]) and (|10p , the 
amplitude of the complex envelope function E(x, t) reads, 



\E(x,t)\ = l^t- + A 2 tanh 2 {p x + 8bpk(k 2 +p 2 )t) 



1/2 

- — h A 2 tanh 2 (p x + 8bpk(k 2 +p 2 )t)\ 

V 

which may approach nonzero when the time variable approaches infinity. 

3. A series of exact solutions for Eq.fll]) by using F-expansion 
method 



We suppose that the solution of ([T|) is of the form 

q{x,t) = P(r)e ij? , r = B(x-ut), r/ = (-kx + ut + 0), 
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(11) 



where P(t) is a real function, and B, u, k, to are real constants to be determined. Substituting 
Eq. ljlip to Eq.([T]) and separating the real and imaginary parts, one may obtain the following 
equations, 

- B{v + 2uk + Abu 3 ) P' + Abn B 3 P'" = 0, (12) 

(u + aK 2 + bn A ) P - cP 3 - cdP 5 - B 2 (a + 66k 2 ) P" + bB A P"" = 0. (13) 

The linear ordinary differential equation (|12p has no solitary wave solutions, thus we have to 
take k = v = 0. In this case Eq. (fT2j) is satisfied identically, and Eq. ([T3|) becomes, 

coP - cP 3 - cdP 5 - aB 2 P" + bB A P"" = 0. (14) 

Eq. (|14p can be solved by using the F-expansion method [3 - [6j. According to the F- 
expansion method, we suppose, 

n 

P(t) =Y J A i F i (r), A n ^0, (15) 

i=0 

where Ai(i = 0, ■ ■ ■ ,n) are real constants to be determined, the integer n is determined by 
balancing the linear highest order term and nonlinear term. And F(t) in (I15D satisfies, 



dF(r) 



2 



ho + h 2 F( T y + h A F{r)\ (16) 



dr 

where ho, h 2 , /14 are real constants. By balancing the linear highest order derivative term P"" 
with nonlinear term P 5 in Eq. (jl4|) . we find n = 1. Thus Eq, (ll5p becomes, 



P(t) = A + A 1 F(t). (17) 

Substituting Eq. (|17p into Eqs. (|14p along with Eq. (|16|) . collecting all terms with the same 
power of F 3 {r){j = 0, • • • ,5), and equating the coefficients of these terms yields a set of 
algebraic equations with respect to Aq-, A\, B, u, a, b, c, d, ho, h 2 , h^. 

dcA^Ao = 0, co A - cA 3 - dcA 5 = 0, 

10<icyli 2 A 3 + 3c^i 2 ^o = 0, 2AbA l B 4 h 4 2 - dcA^ = 0, 

20bA 1 B 4 h 2 h 4 - 10 dcA^Ao 2 - cA x 3 -2A 1 B 2 ah i = 0, 

u Ai - 3 cAiAo 2 - A 1 B 2 ah 2 + 12 bA x B A h A ho + bA x B 4 h 2 2 - 5 dcA x A A = 0. 
Solving the above algebraic equations, we have a set of nontrivial solution, 



A - n A - + !2bB 2 h 4 
A ° - °' Al ~ ± Vd(106S%-a)' 

oj = B 2 (ah 2 - 12bB 2 h 4 h - bB 2 h 2 2 ), c = rf ( 10 ^2-a) 2 



(18) 



Special analytical solutions to Eq. (|16p exists for certain choices of the constants ho, h 2 
and /14. When /i = 1, h 2 = —(1 + m 2 ), = in 2 , Eq. (fT6|) has the solution F(t) = sn(T, m). 
From Eq. ljlip and Eq. (|17p . Eq.([T]) has the Jacobian elliptic sine function solution, 



12B 2 bm 2 , ,■_/>-( n- 1,1 n-i,,,, i 



sn(Bx in) e* ^ B (B b+B bm +14B bm +am + a ) t+e ) 



~\ d{WB 2 b + WB 2 bm 2 + 

where B is determined by d(WbB 2 + Wm 2 bB 2 + a) 2 - 66c = 0. 

When ho = 1 — m 2 , /12 = 2m 2 — 1, /14 = —m 2 , Eq. (jl6p has the solution -F(r) = cn(r, m). 
Inserting it into (I17p and using the transformation ()lip . Eq.([T]) has the Jacobian elliptic 
cosine function solution, 



\W-bl,r cniBx m) <J(B 2 (l<5B 2 bm 2 ~16B 2 b m 4 -B 2 b-a+2am 2 )t+e) 



^ \ d(10 B 2 b - 20 B 2 bm 2 + a) 

where B is determined by d(a + WbB 2 — 20bB 2 m 2 ) 2 — 66c = 0. 

Some solitary wave solutions can be obtained if the modulus in approaches to 1. For 
example, when m — > 1, the solution q 2 {x,t) degenerates to the kink type envelope wave 
solution, 



126 B 2 



t^{Bx)e i{ -- 2B2 ^ B2b+a ^ t+e \ 



\ d(20B 2 b + a) 
where B is determined by d(a + 20bB 2 ) 2 — 66c = 0. 

When m — >• 1, the solution q%{x, t) degenerates to the bell type envelope wave solution, 



q$(x,t) 



V2Dlh -sech^x)^^ 2 ^ 26 )^), 



\| d(a- 10 B 2 b) 

where B is determined by d(a — WbB 2 ) 2 — 66c = 0. 

As pointed out in Ref.[3J, Eq. (|16p has many other Jacobi elliptic function solutions in 
terms of dn(^), ns(£), nd(^), nc(£), sc(£), cs(£), sd(£), ds(^), cd(^), dc(£) as well as the 
corresponding solitary wave and trigonometric function solutions. For simplicity, such types 
of solutions to Eq.([T]) are not listed here. 

With the aid of Maple, we have checked the solutions qj(x,t)(j = 1, ■ ■ ■ ,5) by putting 
them back into Eq.(pQ). 
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